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Abstract. Very few explicit inflationary scenarios are known to generate a large bis- 
pectrum of orthogonal shape. Dirac-Born-Infeld Galileon inflation, in which an induced 
gravity term is added to the DBI action, is one such model. We formulate it in the lan- 
guage of the effective field theory of inflation by identifying the unitary gauge operators 
that govern the behavior of its cosmological fluctuations. We show how to recover rather 
easily from this its power spectrum and bispectrum, which we calculated previously using 
standard cosmological perturbation theory. We push our calculations up to the deter- 
mination of the fourth-order action and of the trispectrum, in which shapes absent in 
k-inflation arise due to the presence of higher-order derivative operators. We finally dis- 
cuss the combined constraints set on this model by current observational bounds on the 
bispectrum and trispectrum. 
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1 Introduction 

The deviation from perfect Gaussian statistics of the primordial fluctuations promises to 
allow a more precise discrimination between competing scenarios of the early Universe 
than has hitherto been possible (see for instance [1, 5] for recent reviews). In this respect, 
one of the important realization of the last decade has been the identification of a dictio- 
nary between, on one side, various mechanisms for generating the primordial curvature 
perturbation £ and, on the other size, various momentum-dependences of its three-point 
correlation function, i.e. various shapes of the primordial bispectrum [6, 7]. The two 
best known shapes, local [8] and equilateral [9], have a well understood physical origin. 
The former is generated classically by the non-linear evolution of the primordial fluctu- 
ations on super-Hubble scales and it is associated to the presence of several degrees of 
freedom during inflation [10, 11] (or other competing early-universe scenarios). The latter 
emerges as the result of non-trivial quantum correlations between fluctuations that cross 
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the Hubble radius at about the same time during inflation, and it is generated in models 
with non-standard kinetic terms, single- or multi- field [12-15]. Archetypal scenarios that 
represent these two classes of models are hence known, like the curvaton [16] and DBI 
inflation [12, 17] respectively 1 . Yet another type of non-Gaussianities, namely orthogonal 
ones [19], is often considered in the literature (see for instance the recent constraints from 
WMAP9 [20]). However, its status is rather different from the two others described above. 
It has been constructed in the formalism of the effective field theory of inflation [21] as 
the outcome of possible partial cancellations between equilateral-type shapes generated in 
single-field models, without any motivation from a particular model, but rather with data 
analysis in mind. 

While this is a perfectly valid origin, we find it important to have at hand concrete mi- 
crophysical early-universe scenarios able to generate a large orthogonal bispectrum, all the 
more so as the WMAP nine-year analysis indicated a (non-significant) 2. 45a hint of such 
a signal [20]. Dirac-Born-Infeld Galileon (DBI Galileon) inflation, considered in references 
[22, 23], provides the first example of such a model 2 . Its origin is physically transparent 
and its action can be written in simple geometrical terms: it consists in supplementing the 
DBI action, describing the motion of a probe D3-brane, with an induced gravity term. It 
can hence effectively be seen as an extension of the well known DBI inflationary scenario 
with only one new dimensionless parameter setting the strength of the induced gravity. 
When varying the latter, the shape of the bispectrum varies continuously, interpolating 
between highly correlated and highly anti-correlated with the equilateral template, while 
it is of the orthogonal type in a transitionary region. 

The purpose of the present paper is two-fold. As the orthogonal shape has been first 
motivated and constructed in the effective field theory formalism, we would like to formu- 
late the DBI Galileon model in this language, i.e. to identify the operators and related 
mass scales entering into the Lagrangian describing its fluctuations. As we will see, the 
way the orthogonal shape appears in this model is similar to, but different in details from, 
the one in its original construction in [19]. Our second aim is to investigate the possibility 
of observationally falsifying DBI Galileon inflation as an inflationary scenario candidate 
for explaining the large orthogonal bispectrum suggested in WMAP data [20] . Indeed, we 
recently pointed out the fact that the trispectrum can be used, already with current data, 
to constrain such candidates in the simplest theoretical framework [25] . Although the two 
setups are different, we will see that similar conclusions can also be reached here, and that 

1 Naturally, these two shapes do not constitute separate possibilities and a linear combination of them 
can be present in the same model, like in multifield DBI inflation [18]. 

2 Interestingly, a bispectrum with a significant overlap with the orthogonal shape was recently shown 
to arise naturally in a different setup [24] . 
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the trispectrum is putting the DBI Galileon model under strong observational pressure, 
at least in the la region favored by WMAP9 constraints on the bispectrum. 

The plan of our paper is the following: in section 2, we introduce the DBI Galileon 
model and briefly discuss its background evolution. In section 3, we formulate it in the 
language of the effective field theory of inflation. We show in particular how to use 
the latter to recover rather easily the leading-order quadratic and cubic actions in the 
cosmological fluctuations, and hence the power spectrum an bispectrum determined in 
[22, 23]. We also derive the action quadratic in fluctuations, and the primordial trispectrum 
generated in this model. The section 4 is then devoted to a discussion of its observational 
status. We conclude in section 5 and collect technical details, lengthy results and some 
plots in three appendices. 

2 DBI Galileon: action and background 

In this section, we briefly review the aspects of the DBI Galileon model that are relevant 
for what follows, referring the reader to [22, 23] for an extensive study of it, in particular 
its multifield aspects. For simplicity, we consider only its single-field version here, or 
equivalently we neglect the influence of the entropic perturbations on the adiabatic one. 
The action we consider reads 



where the first term is the standard Einstein-Hilbert action involving the cosmological 
metric g^ u , the second term is the Dirac-Born-Infeld action [17] 



and the third term is proportional to the Einstein-Hilbert Lagrangian of the so-called 
induced metric 3 



We also consider a constant warp-factor / for simplicity (see [23] for a discussion of the 
more general case), acknowledging that the effects of a non-constant / would anyway be 
negligible in a slow-varying approximation. 

3 In the higher-dimensional framework which gives rises to the effective four-dimensional Lagrangian 
(2.1), the metric (2.3) is proportional to the induced metric on the brane, hence its name. 
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9nv = 9tiu + f dn<f>d v <f> . 
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In a spatially flat Friedmann-Lemaitre- Robertson- Walker (FLRW) background space- 
time, of metric 

g^dx^dx" = -dt 2 + a 2 (t)dx 2 , (2.4) 

the background equations of motion derived from the action (2.1) can be cast in the simple 
form 4 [22, 23] 

atfta + D^ + I^-i) (2.5) 



where e = —H/H 2 denotes the standard inflationary parameter, with H = a/a being the 
Hubble parameter, 

4 = l-/0 2 , (2.7) 

and where s = c-t>/{Hc-p) 5 ■ Like in standard brane inflation, the value of cfj distinguishes 
the slow-roll regime, where c|> ~ 1 and the brane action (2.2) takes a canonical form, from 
the so-called relativistic, or DBI, regime, where the inflaton almost saturates its speed 
limit (<fi 2 < 1/f) and c|, <C 1. It is this regime, where the non-linear structure of the 
action (2.2) must be taken into account and large non-Gaussianities are generated, that 
we are particularly interested in. From Eqs. (2.5)-(2.6), one can show that achieving a 
phase of quasi de-Sitter expansion - e <C 1 - in the relativistic regime, requires, in addition 
to the usual condition to enter into the DBI regime 

cvfV > 1 , (2.8) 

that 

M 2 < c%Ml x . (2.9) 

Despite this restrictive condition on the mass scale associated to the induced gravity, the 
effects of the latter on the cosmological fluctuations can actually be significant, and are 
conveniently measured by the dimensionless parameter 

fH 2 M 2 

a = - — 2 . (2.10) 



4 As usual in single-field inflation, the equation of motion derived from varying the action (2.1) with 
respect to cj) is a consequence of the gravitational equations (2.5)-(2.6). 

5 We use the symbol cf,, and not the more conventional c 2 s , because, contrary to standard brane inflation, 
cx> does not coincide here with the speed of sound of scalar perturbations. 
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The latter is not necessarily small during inflation, although we will see it has to be less 
than 1/9 to avoid the presence of ghosts. In the inflationary regime of interest where the 
conditions (2.8) and (2.9) are fulfilled, note eventually that Eq. (2.6) reduces to 



[ i>i 



At leading order in a slow- varying approximation, we can consider cv and a as constant, 
and we will see in section 4 how observational constraints on the primordial bispectrum 
and trispectrum directly translate into constraints on these two parameters. 



3 DBI Galileon: fluctuations in the effective field theory 

In references [22, 23], we calculated the second- and third-order action of cosmo logical 
fluctuations in DBI Galileon inflation, from which we deduced the corresponding power 
spectrum and bispectrum of the curvature perturbation (. In this section, we show how 
these results can be recovered in an elegant way in the effective field theory formalism. 
Moreover, the latter is powerful enough that we are able to derive rather easily the fourth- 
order action, as necessary for the study of the primordial trispectrum. 

We refer the reader to [21, 26] for details about the construction of the effective field 
theory of inflation, only summarizing here the strategy we followed in our calculation: we 
first expressed the action (2.1) in the unitary gauge in which the scalar field <fi is unper- 
turbed (subsection 3.1), paying a particular attention to put our result in a convenient 
form for the second part of the calculation (subsection 3.2), where we used the Stiickelberg 
trick to explicitly reintroduce the scalar perturbation tc and derive the action governing 
its dynamics in the decoupling regime in which the mixing with gravity is irrelevant. 

3.1 Action in the unitary gauge 

We set ourselves in the unitary gauge in which the scalar field is unperturbed, i.e. such 
that <p = <f)(t), and use the ADM formalism [27] in which the cosmological metric is written 
in the form 

g^dx^dx" = -N 2 dt 2 + h^dx 1 + N i dt)(dx j + N j dt) , (3.1) 

where hij is the spatial metric, N l is the shift vector and N the lapse function. In the 
unitary gauge, the induced metric (2.3) takes the very simple form 

g ilv dx ll dx v = -N 2 dt 2 + h ij (dx i + N l dt)(dx j + N j dt) (3.2) 

where only the lapse function 

N 2 = N 2 - 1 + 4 (3.3) 



- 5 - 



is modified compared to Eq. (3.1). Using the ADM form of the Einstein-Hilbert action, 
the induced gravity term S- m( i = / d 4 Xy/^gR[g] reads 

S M = ^-Jd*x^g (^RW + | ( Kij K- - K 2 ) j (3.4) 

where R^ is the scalar Ricci curvature of the spatial metric and Ky = ^ [hij — 2V(jiVj) 
is the extrinsic curvature of constant time hypersurfaces. In the effective field theory of 
inflation as formulated in [21], the 3d scalar Ricci curvature R^ is traded for the four- 
dimensional one R using Gauss-Codazzi relations. While this is useful for classification 
purposes in a general formalism, it turns out here to be easier to use the two variables for 
practical calculations, as also noted in the different context of the effective field theory of 
dark energy [28]. Using [21] 

R {3) =R + K 2 - K, tJ K i] - 2V '^V ' v n v ) + 2V I/ (n"V (1 n i/ ) , (3.5) 

where is the unit vector perpendicular to constant time hypersurfaces (we choose it 
such that n° > 0), and the fact that g 00 = —1/N 2 , we therefore conveniently write 

/ ^x^g^R^ = J d'x^-g (j\ + (1 - c 2 v )g™ - c^j R^ 

+ J d 4 x^g~c v (R + K 2 - KijK^ - 2V M KV,rO) (3.6) 
where in the last line we used the fact that [21] 

d^v^M^V.KVX) = - / d 4 xy/=gd v (cj>(t)) b"V/ = (3.7) 



as d v (cx>(t)) oc n u . While the advantage of the mixed form Eq. (3.6) may not be obvious 
at this stage, it will become transparent in the next subsection where we will explicitly 
reintroduce the scalar fluctuation n. Using Eqs. (3.4) and (3.6), we can finally rewrite the 
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full action (2.1) in the form (convenient for what follows) 



S 



+ 



+ 



d 4 x 

M 2 

~Y 

M 2 



cv 



M 2 
M 2 , 



R 



d 4 x 



+ M 



+ (1-4)^°° 
d 4 Xv^ (j/l + (1 - 4)0«> _ c ^ ^(3) 

d 4 Xv^n^(c c (t))V^ 



c c (^-^ - if 2 ) 



+ / &Xy/=g 



1 

7 



4)9 



00 



(3.J 



3.2 Reintroducing the Goldstone boson 



We aim at deriving the action governing the dynamics of the cosmological fluctuations 
about the background discussed in section 2. For this purpose, and following [21, 28, 29], 
we separate Eq. (3.8) into 



S 



M 2 , 



g(t)R-A(t)-c(t)g 



oo 



_|_ g(n>2) 



(3.9) 



where S^™- 2 ) explicitly starts quadratic in the fluctuations Sg 00 = g 00 +l, 5Kij = Kij—Hhij 
and SR^ = R {3) . As explained in [21], g(t) can be set to one by a conformal transformation 
of g^ v . While we could follow this strategy here, we find it more convenient to keep a non- 
trivial g(t), as arising naturally from the first line in Eq. (3.8). In this sense, although we 
are dealing with a phase of single-clock inflation, we are actually using the tools of the 
effective field theory of dark energy developed in [28, 29]. We leave the explicit form of 
git), A(i) and c(t) to the appendix A, where we also give details of the intermediate steps 
that lead to the result S {n ^ 2) = J d A Xyf^gC {n ^ 2) with 



£(n>2) 

(1 - 15a 
16/cf 



l-9a) 



(c^ 00 ) 2 + ^-8g m 8K + ^-{bKifiK 13 - 5K 



HM 2 



,00 



M 2 



8/c 3 



M 2 



v 



'D 



.^y- 3I ^l(5g 00 ) 2 5K- 



2c v 
M 2 



Sg^iSKijdK 13 - 5K 2 ) - 



'D 

5(1 - 21a) 
128/4 



(Sg 



00\4 



+ 



4c|, 
5HM 2 



4c|, 
3M 



4c© 
M 2 



8g 00 8RW (3.10) 



164 



(5g°°) 2 5R^ 
M 2 



(Sg 00 fSK + ^(Sg 00 ) 2 (SK^ - 5K 2 ) + ^fo^SR™ 



164 



32c^ 



up to quartic order in fluctuations. Here, we kept leading order terms in the relativistic 
regime 4^1 an d i n a slow-varying approximation. In particular, we neglected terms 
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involving time-derivatives of (it is easy to take them into account if needed), arising for 
instance from the fourth line in Eq. (3.8). This is an important result of our paper: the 
identification of the operators and related mass scales that specify the DBI Galileon model 
in the effective field theory of inflation in the unitary gauge. The first line (respectively 
second and third) contains the terms that start explicitly quadratic (respectively cubic 
and quartic) in the fluctuations. Note that the result simplifies considerably for the DBI 
model, corresponding to M 2 = a = 0, and for which the only present operators are powers 
of 5g 00 . 

While Eq. (3.10) is an important step in our calculation, the true effectiveness of 
the effective field theory of inflation, however, relies on the gravitational analogue of the 
equivalence theorem for the longitudinal components of a massive gauge boson [30]: upon 
performing a time diffeomorphism t — > t + ir(x) on the gauge-fixed action (3.8)-(3.9), thus 
restoring full time-diffeomorphism invariance 6 , one can show that the scalar fluctuation 
7r reintroduced this way decouples from the gravitational sector at high enough energies, 
allowing to neglect the complications of the mixing with gravity. This latter fact was also 
shown to be valid in the relativistic regime of DBI Galileon inflation in [22, 23], and we 
consider only this decoupling regime in what follows. 

The effect of the Stiickelberg diffeomorphism t — > t + tt(x) on explicit functions of 
time is h(t) h(t + it) « h(t) + h(t)n + ^vr 2 + . . . The terms in 7r generated in this 
process are negligible in the slow-varying approximation and we don't keep track of them, 
consistently with the fact that we discarded terms involving derivatives of c*r>- Restoring 
7T thus only affects the operators through the replacements [21, 28, 29] 

6g oo _ 2 7T - tt 2 + ^ , (3.11) 

cr 

SKij — > —Hnhij — didjix , (3-12) 

5K -»■ -3Htt J , (3.13) 

cr 

5R {3) ^0, (3.14) 

again, in the decoupling regime in which the metric is left unperturbed. Note eventually 
that 7r is related to the curvature perturbation ( by the simple relation ( = —Hit at linear 
order and at leading order in a slow- varying approximation [21]. 

3.2.1 Quadratic action and power spectrum 

From the above formulae, it is straightforward to derive the main second-order action 



5( 2 ) = / dtd 3 xa 3 



1 — 9a) . 2 (1 — 5a) (<9j7r 



,2 



-7T 



2fc 3 v 2fc v a 2 



(3.15) 



6 Thc field 7r is assumed to transform under the time diffeomorphism i — t + ^(x^) as tt = tt — ^°(x M ), 
such that t + tt is invariant. 
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in agreement with the one found in [22, 23], and where we used the integration by part 
J dtd 3 x a 3 — J dtd 3 x a y(<9j7r) 2 . As mentioned in section 2, tt becomes a ghost when 
a > 1/9, hence our limitation to the range < a < 1/9. From Eq. (3.15), one deduces 
the speed of sound c s at which the curvature perturbation propagates: 



1 — 5a 

"•^Vi^' (3 ' 16) 

Note that for any cx>, there is a limiting value of a above which c s in Eq. (3.16) is greater 
than one. However, in the relativistic regime in which our calculations are applicable, for 
od < 0.1, these values are very close to 1/9. We will see that the theory becomes strongly 
coupled as a 4 1/9, not surprisingly as the kinetic term in (3.15) goes to zero in this 
limit. For practical purposes, there is thus no issue of superluminality, as the region where 
it arises is already excluded by observational bounds on non-Gaussianities. 

In the almost de-Sitter regime of interest, in which a ~ —1/(Ht) with r the conformal 
time such that dt = adr, the linear equation of motion derived from Eq. (3.15) is simply 

n + 3Hn-c 2 s — £ =0, (3.17) 
a z 

whose standard solution with Bunch-Davies normalization and initial conditions reads (in 
Fourier space) 

7T k (r) = J ]~^* _L — (1 + ikc s r)e~^ T . (3.18) 
V 1 - 5a y/4ec s k 3 Mpi 

We then easily deduce the power spectrum of the curvature perturbation: 

k 3 r> , v„ 1 - 3a 1 H 2 . 
"C- ^I6(^0)|' = — (3.19) 

3.2.2 Cubic action and bispectrum 

Around the time of Hubble crossing kc s r ~ —1, it is straightforward to verify from 
Eq. (3.18) that 7T ~ Hit and d^ja ~ Htc/c-d. Taking this into account (note in par- 
ticular that amongst terms with the same number of derivatives, the ones with the highest 
number of space derivatives dominate in the relativistic regime c|, <C 1), we deduce from 
Eqs. (3.10) and (3.11)-(3.14) the leading-order cubic action: 



S( 3 ) = / dtd 3 xa 3 



l-15a). 3 C 1 - 9 «) • (W J 

-7T r-5 TT- 



2/4 2/4 a 2 

3HM 2 , 2 d 2 n 3HM 2 {d t n) 2 d 2 7r 



-TT 

4 a 2 44 a 4 



(3.20) 
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where we made use of 



/ 



dtd 3 x - ((didjTr) 2 - (<9 2 tt) 2 ) = [ dtd^x^d^Yd 2 ^ , (3.21) 



derived by successive simple integrations by part. Note that the cubic action (3.20) con- 
tains terms involving second-order (spatial) derivatives of it (in its second line) as a result 
of the presence of extrinsic curvature terms in the unitary gauge action (3.10). However, 
as it was explained in a general context in references [31-33], it is legitimate to simplify 
the interacting action of cosmological fluctuations by using the linear equation of motion 
derived from the quadratic action. Using this fact, one can actually trade the two higher 
derivative operators 7r 2 <9 2 7r and (<9;7r) 2 <9 2 7r for the single-derivative ones 7r 3 and 7r(<9j7r) 2 (see 
for instance [34] and Eqs. 4.5 and 4.6 in [23]). The effective third-order action then reads 



= I dtd A xa 



cS - , ..... 



v 



•3 2 (di^) 

A^sTT - A^Q^CpTT y 

(X 



(3.22) 



with 



A*3 = 1 - 3a(5 - 4A 2 + A 4 ) , (3.23) 
A Hd ^)2 = 1 - a(9 - 3A 2 ) , (3.24) 

c s \ 1 — ba 

We thus recovered the cubic action derived in [22, 23], albeit with considerably less effort 
thanks to the effective field theory formalism. The primordial bispectrum deduced from 
Eq. (3.22) reads [22, 23] 



3 



(C(fcx)C(fc 2 )C(fc3)> = (2k?6{£ ^V f^^ (3-26) 

where is the primordial power spectrum given in Eq. (3.19) and the shape function 

11/ A 2 A 

S(ki, k 2 , h) = ~2^2" i - Q a ySA^S^ + —A^g.^S^g.^j (3.27) 

is a linear combination (function of a) of 

= ^ (3.28) 

S m7T)2 = ~^jr^ k i ■ fc 2 (2hk 2 - k 3 K + 2K 2 ) + 2 perms. (3.29) 

where K = k\ + k 2 + k$. The bispectrum generated in DBI Galileon inflation can not 
be mapped to the one in the simplest effective field theory of inflation [21]. Yet, its 



- 10 - 



qualitative features are similar: its overall amplitude is fixed by the speed of sound c s (or 
equivalently by cp) while its shape depends continuously on a dimensionless parameter 
(here a). In particular, the former interpolates between highly anti-correlated (respectively 
correlated) with the equilateral template as a — > (respectively a — > 1/9), while being of 
the orthogonal type in a transitionary region centered around a = 0.097. This is shown 
in Fig. 1 where we plot the scalar product, in the sense defined in [7], between the DBI 
Galileon shape (3.27) and the equilateral and orthogonal shapes as a function of a (see 
[23] for more details and section 4 for observational constraints). 

1.0 I ^ 1 




0.105 



Figure 1. Correlation of (minus) the DBI Galileon shape (3.27) with the equilateral template 
(blue) and the orthogonal template (red) as a function of a. 



3.2.3 Quartic action and trispectrum 

Following the same logic that led to the cubic action (3.20) from Eqs. (3.10) and (3.11)- 
(3.14), we find the leading-order quartic action 



S (4) 



dt d x a 



5(1 -21a). 4 3(1 - 15a) . 2 (^tt) 2 (1 - 9a) ((cV 

'7T r^r^ VT h 



,2\2 



8/4 



v 



4/c 5 



v 



a/ 



8/c 3 



v 



a 1 



+ 



M 2 
44 



o • 2 2 

Sir — c£,- 



(cVr) 2 \ {d l d j n) 2 - (cfTr) 2 5HM 2 . 3 «9 2 vr 3HM 2 . (<9 j7 r) 2 <9 2 7r 



7T 



-7T- 



'D 



'D 



(3.30) 



Similarly to the case of the cubic action, one can simplify (3.30) by using the linear equation 
of motion (3.17), leading to (see appendix B for intermediate steps): 



S. 



(4) 
cff 



dt d 3 x a 



8/4, 



5 j I — n ( 21 - 18A 2 + yA 4 



7T 4 



l-9a)c v — 



,2\2 



6 (1 - 5a (3 - 2A 2 )) c^^f- + ^ ( 3 vr 2 - " ^ 



a 4 



(3.31) 



which can not be further simplified, i.e. where the quartic operators are independent of 
each other. This fourth-order action can be usefully compared to its counterpart in k- 
inflation [35] (see also [36]). One observes that the two operators in n 2 ((didjir) 2 — (dfir) 2 ) 
and (dkTr) 2 ((did j-n) 2 — (dfir) 2 ), which, as we stress again, can not be traded for the three 
other ones, are not generated in k-inflation. This is to be contrasted with the cubic action 
(3.22) and is a direct consequence of the presence of higher-derivative operators involving 
extrinsic curvature terms in the unitary gauge action (3.10). 

We find then that the primordial trispectrum generated in DBI Galileon inflation 
reads 



(c(fc 1 )c(fc 2 )c(fc3)c(fc 4 )) c = (2^5^ fe ^ c 3 r(fci, ^fc^3 i2,fci4) ( 3 - 32 ) 

i=l 

with kij = \ki + kj\ and 

r = \ ( J2 A a T ci + Yl ( 3 - 33 ) 

C V \i=l j= l J 

being the sum of 8 contributions: 5 contact-interaction trispectra T C j - corresponding to 
the 5 quartic operators in (3.31) - and 3 scalar-exchange trispectra T s j, solely determined 
by the cubic (and quadratic) action (3.22). We use the same notations as in [35] for 
T c i,T c2 ,T c3 ,T sl ,T s2 and T s3 , whose explicit and lengthy expressions can be found there, 
while we defined 

T C 4 = 



T, 




+ 23 perms. (3.34) 

, 5 = (k 1 -k 2 )((k 3 -k 4 ) 2 -A;3 2 A: 4 2 ) x 

%^ + 90%^ ) + 23 perms. (3.35) 
K 6 iv* J 

for the trispectra generated respectively by the higher-derivative operators ^((didjix) 2 — 
(dfir) 2 ) and (dk7r) 2 ((didjn) 2 — (df-ir) 2 ) (here K = k% + k 2 + k 3 + k A and the 'perms.' stands 
for the permutations of the 4 wavevectors kj). We leave the explicit expressions of the 
dimensionless coefficients A ci and A S j (functions of a) to the appendix B, where we also 
give details about the calculations leading to the result (3.33). Besides the appearance of 
the new shapes (3.34) and (3.35), let us stress one important difference compared to the 
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trispectrum generated in k-inflation. There, the specific combination of the coefficients of 
the cubic action and of ((<9j7r) 2 ) 2 and n 2 (diii) 2 are such that these two operators cancel at 
leading order in the quartic Hamiltonian. Such cancellations do not arise in our model, 
where T C 2 and T C 3 contribute at leading-order (i.e. as l/cf>) to the trispectrum (3.33). 

4 Observational constraints 

In the simplest setup of the effective field theory of inflation, computationally equivalent 
to k-inflation, the orthogonal shape of the bispectrum arises due to partial cancellations 
between equilateral-type shapes [19]. This implies that generating an orthogonal bispec- 
trum at the level suggested in WMAP data requires in this context a speed of sound of 
inflaton perturbations c s as low as about 0.01. Unless inflation occurred in the region of 
parameter space where similar partial cancellations that leads to the orthogonal bispec- 
trum arises at the level of the trispectrum, the amplitude of the latter, scaling as 1/cf, 
is of order 10 8 and hence is comparable to (and greater than) current observational con- 
straint f£ L = (-3.11 ± 7.5) x 10 6 [37]. DBI Galileon inflation does not belong to the 
simple of class of models in which we quantitatively developed this argument in [25] , and 
the results derived there can not be simply applied here. Yet, we have seen that the two 
frameworks share similar qualitative features and it is therefore reasonable to expect sim- 
ilar constraints on this model from the trispectrum. We investigate this question in what 
follows, equipped with the predictions of the preceding section. 

4.1 Constraints from the primordial bispectrum 

In [20] , the WMAP team translated observational constraints on the primordial bispectrum 
to constraints on the two parameters specifying the cubic action of the simplest effective 
field theory of inflation. As the bispectrum (3.27) generated in DBI Galileon inflation is 
similar to the one generated in this framework (it is a linear combination of the same two 
shapes), it is trivial to perform the same analysis here. In particular, Eq. (57) in [20] 
directly translates here to 

f£ L = 2n 1 n , (-0.276 A 2 - 0.0785 A,*) (4.1) 

C-p \ 1 — i)Q ) 

f°NL = ^ \ Qa) (0.0157 A 2 Av ( ,9 t7r )2 + 0.0163 A^) , (4.2) 

where A^, A^g.^2 and A are given in Eqs. (3.23)-(3.25). This provides the link between the 
two dimensionless parameters (cd, a) of our model and the two parameters (f^ L , /atl 1 ) as 
constrained by WMAP 9. Following exactly the same logic as the one in [20], we then easily 
deduce the la, 2a and 3a confidence regions displayed in Fig. 2. They are qualitatively 
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Figure 2. la (red), 2a (orange) and 3a (blue) confidence regions in the (ep,a) plane derived 
from the WMAP9 bispectrum constraints [20]. 

similar to the constraints on the parameters of the simplest effective field theory of inflation 
(c s ,A) given in [20]. In particular, the la region favored by data (in red) corresponds to 
low values of of about 0.01 and to the narrow region 0.096 < a < 0.098 where the 
bispectrum (3.27) significantly overlaps with the orthogonal template (see Fig. 1). The 
fact that a negative is favored at 2a translates at this level of confidence into the 
lower bound a > 0.095 (in orange). Note also that as a — > 1/9, the amplitude of fff L and 
/jv£ h in Eqs. (4.1)-(4.2) grow unboundedly, corresponding to a strongly coupled region. 
This explains the sharpness of the bound on cx> in the upper half region. Conversely, for 
a < 0.090, the bispectrum (3.27) is of equilateral type and is nearly insensitive to the 
strength of the induced gravity. 

4.2 Constraints from the primordial trispectrum 

A simple measure of the amplitude of the trispectrum is given by the parameter t^i 
defined such that [35] 

— T(h, k 2 , k 3 , fc 4 , fcu, k u ) > t NL , (4.3) 

K 6 limit 
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where RT stands for the regular tetrahedron limit k% = k 2 = k 3 = k^ = fc 12 = A; 14 = k. 
While the evaluation at this single point might seem rather ad-hoc, it can be motivated by 
the fact that trispectra of quantum origin generated around the time of Hubble crossing are 
generally the largest for the configurations where both the external and internal momenta 
are of similar magnitude, i. e. near this regular tetrahedron limit. This was shown explicitly 
in reference [35] for the trispectra T c i,T c2 ,T c3 ,T s i,T s2 and T s3 , and this can easily be 
checked to hold as well for T c4 and T c5 in Eqs. (3.34)-(3.35). 



Using this similarity, the observational constraint relevant for our purpose, t 



cq 
NL 



(—3.11 ± 7.5) x 10 6 , was derived in [37] (based on methods developed in [38]) under the 
assumption that the trispectrum shape T in Eq. (3.32) can be well approximated by the 
simple and representative 'equilateral-type' shape T cl generated by the operator 7r 4 : 



{k\k 2 k 3 k, 



\2 



T ci = 36 ,, v , J , 7 , ._ . (4.4) 
(fci + k 2 + k 3 + k A f y ' 

More precisely, this constraint - the only one to date on a trispectrum from quantum 
origin - puts a bound on the parameter f^ L defined such that T = t e ^ L /tNL{T c i) x T c \. It 
is therefore only applicable in our model for values of a such that the trispectrum (3.33) 
can be well approximated by T c \. Fortunately, we have checked this to be the case in 
the interesting region a > 0.095 favored (at 2a) by the observational constraints on the 
bispectrum (see Fig. 2 in the preceding subsection) 7 . For this purpose, we used the same 
strategy as in [25] and represented 2-dimensional slices of our trispectrum (3.33) (as a 
function of a) in different representative limits, as was originally done in [35]. We display 
the outcome of this procedure for the representative value a = 0.097 (in the middle of 
the la region favored by the bispectrum) in appendix C, from which the approximate 
proportionality with T cl is manifest. 

To summarize, the possibility - highlighted in a different context in [25] - of par- 
tial cancellations leading to the emergence of an 'orthogonal-type' trispectrum qualita- 
tively different from T cl does not arise in DBI Galileon inflation in the interesting range 
a > 0.095. The estimator (4.3) can hence be reliably compared in this regime to the ob- 
servational bound fj^ L = (—3.11 ± 7.5) x 10 6 (68 % CL). We show in Fig. 3 the confidence 
level at which DBI Galileon is excluded (or not) by this trispectrum constraint, i.e. the 
represented quantity is \tNL — t% L \/ At^L, where t* NL = —3.11 x 10 6 and AtNL = 7.5 x 10 6 . 
The upper figure (respectively the lower figure) displays the whole 3a region (respectively 
zooms in the la region) favored by the bispectrum constraints, for a > 0.095. The conclu- 
sion is clear: within the la region favored by observational constraints on the bispectrum, 
most of the parameter space (c©, a) is already excluded by the observational bound on 



7 Note also that for a < 0.095, allowed values of Cx> from bispectrum constraints are such that the 
amplitude of the trispectrum is too low to be constrained by data. 
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the trispectrum, while the amplitude of the latter is too small outside this region to be 
efficiently constrained by data. 

5 Conclusions 

The DBI Galileon model, in which an induced gravity term is added to the DBI action, 
is one of the few known explicit inflationary scenarios able to generate a large primordial 
bispectrum with a significant overlap with the orthogonal template. As the latter was 
first motivated in the formalism of the effective field theory of inflation, we formulated 
DBI Galileon in this language. In particular, we identified the various operators and re- 
lated mass scales that specify the behavior of its fluctuations in the unitary gauge. As a 
consequence of the presence of higher-derivative operators, we have seen that the way the 
orthogonal shape arises here differs in details from the one in its original construction. Fur- 
thermore, we determined the leading-order action quartic in fluctuations, and calculated 
the corresponding trispectrum. In the region of parameter space in which the orthogonal- 
type bispectrum is generated, we found that the former can be very well approximated 
by the simple 'equilateral' trispectrum used in current data analysis. Translating the 
corresponding observational bound on constraints on the parameters of the model, we 
showed the combined consideration of the bi- and trispectrum to be very discriminating. 
In particular, it is only in a small window of parameter space that DBI Galileon inflation 
can generate a bispectrum within the la region favored by WMAP9 data without at the 
same time generating too large a trispectrum in conflict with observational bounds. It will 
certainly be interesting to compare the predictions derived here to the results of future 
observations, and in particular to the forthcoming ones of Planck. 

On a more theoretical side, we believe several aspects our calculations to be interesting 
in their own right. In particular, although we considered an inflationary model, we actually 
made use of tools developed in the context of the effective field theory of dark energy 
[28, 29]. The reason behind this is that some aspects of the effective field theory of inflation, 
as formulated in [21], are useful for classification purposes in a general formalism, but may 
less so for practical calculations in concrete models. For instance, although it is always 
possible, through a conformal transformation of the metric, to set to one the function in 
front of the Ricci scalar in Eq. (3.9), it turned out to be simpler not to do so here, related to 
the fact that there is no Einstein frame in the DBI Galileon model. Similarly, although one 
can always trade the 3d Riemann tensor and covariant derivative for their 4d counterparts, 
keeping both type of quantities greatly simplified our calculations. Eventually, note that, 
by using the linear equation of motion to simplify the interacting action, we showed that 
the presence of higher-derivative operators in the unitary gauge action Eq. (3.10) manifests 
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itself as higher-derivative ir interactions at the level of the quartic action only, and it would 
be interesting to study this in more generality. 
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A Calculations in the effective field theory formalism 

In this appendix, we give some details about the calculations performed in subsection 3.2. 



The derivation of Eqs. (3.9) and (3.10) from the action Eq. (3.8) follows the steps described 
in [21] to show that any single-clock model of inflation can be cast in the form (3.9). In 
more details: the way we wrote down Eq. (3.8) is such that its first line is already in the 
desired form. As has no background value, i.e. 8R^ 3 ' = R^ 3 \ the third line is such that 
it only contributes to S^ n - 2 \ as 5R^ times powers of 5g 00 . The last line, corresponding 
to the DBI scenario, without induced gravity, can be simply expanded in powers of 5g 00 . 
For the second line, we conveniently write 



KijK ij - K 2 = 6H 2 - AHK + 6X^5 K ij - (SKf 



(A.l) 



where terms linear in K, here and coming from the fourth line, are dealt with using 



d A x\f—gn^duh 



d^x^g ^-g m h . 



which is in a form ready to be expanded in power of 5g 00 . Following this, we easily find 
Eqs. (3.9) and (3.10) with 



g{t) = l + c v 



p 



(A.2) 



A(t) = V+^(l-c v ) 2 



2fc v 



+ 



M 2 



c(t) 



2 

cv ( 1 



c v + 3Hc v - 2 ( H [ c v 

cv 



3H 2 



;i + 24)(i-4) 



T> 



(A.3) 



2/ VI 



1 + 



M 2 



-c v + 3Hc v + 2 ( H | — - c P 



1 



3H 2 (l-c 



V) 



-v 



(A.4) 
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where we made no approximation. As explained in [29] (see also [28] for an equivalent 
discussion), the gravitational equations of motions derived from an action of the type (3.9) 
reads, in a spatially flat FLRW background: 

c-l&gf-H-ll + ?-*), (A.5) 

T2 Jrj , OZT2 , 1 '9 , $ H 9 



A = M^g{H + 3H 2 + + . (A.6) 

We have verified indeed that these two equations, with g(t), A(t) and c(t) given in 
Eqs. (A. 2), (A. 3) and (A. 4) respectively, are equivalent to the ones given in (2.5)-(2.6). 
This provides a good consistency check of our calculations. 

B Calculations of the quartic action and of the trispectrum 

In this appendix, we give some details about the calculations of the quartic action and of 
the trispectrum. 

• Simplification of the quartic action: 

Between Eq. (3.30) and Eq. (3.31), we traded the operators in n 3 d 2 n and 7r(diii) 2 d 2 iT 
for well known ones generated in k-inflation. For this purpose, we used the (leading- 
order) linear equation of motion (3.17) to replace instances of dfn/a 2 by (a 3 7r)' /(a 3 c 2 ). 
Integrating by parts in time, one then finds 

dtd 3 xa 3 Tr 3 ^= [ dtd 3 xa 3 ^\ (B.l) 
a? J Ac 2 

Using this result, together with 



dtd 3 xa 3 7r 2 ^^ = -- dtd 3 xa 3 7r 3 ^, (B.2) 



, d l irdi7r 1 f , , q o . o d 2 ir 
— -i- = — / dtd 3 xa 3 ir 3 ^- 
a 2 3 J a 2 

obtained by simple spatial integrations by part, one similarly finds 

fi-^JW !* _ /\ + * 3 f^HJd^f , 3H^ 



dtd 3 xa 3 7r- = / dt<Txa 6 — o ^" — ~<r~ + inr* • ( B -3) 

a 4 J \2c 2 s a 2 4c 4 J K ' 

• Calculation of the trispectrum: 

Once having determined the second-, third-, and fourth-order action of cosmological per- 
turbations, here in Eqs. (3.15), (3.22) and (3.31) respectively, the procedure for calculating 
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the primordial trispectrum is well known. It is based on the In-In (also known as Keldysh- 
Schwinger) formalism [39], according to which the expectation value of any operator Q(t) 
in the interacting vacuum of the theory reads 



«?(*)> = (o| 



Texp 



/ / Hj(t')dt' 

— oo(l+ie) 



Q\t) Texp i-i / Hj(lf)df) |0>, 

V J-oo(l-ie) /J 

(B.4) 

where |0) denotes the vacuum of the free theory, T the time-ordering product, Hj the 
interacting Hamiltonian, Q J {t) is in the interacting picture, and the ie indicates that the 
time integration contour should be slightly rotated into the imaginary plane in order to 
project onto the interacting vacuum [10]. While the cubic Hamiltonian is simply 
— L( 3 ), this simple relation does not hold at quartic order. Instead, a generic Lagrangian 
of the form 



£2 = fo^ 2 + 32 

£4 = hoTT 4 + /li7T 3 + h 2 7T 2 + h 3 ( + J 4 



5-07T 3 + gi 7C 2 + g 2 7T + j 3 , 



(B.5) 
(B.6) 
(B.7) 



gives the fourth-order interaction Hamiltonian density [40]: 



si 



where /, g, h and j's are functions of n, diix and £, and the subscripts denote 
of Ti. In our context, this gives 



(B.8) 
the orders 



a 3 eH 2 M^ 1 - 5a 1 



4c? 



+ 



1 c s 1 . _> 

3^2 (^TT, 



3a c|> L 



-4Aci 7r| + 4A c2 cf7r| 



n ^c3 C g -J- 



2\2 



3d A 



(^A (didjTT) 2 - (dfr) 



c5 



21 



(B.9) 
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where 



A c2 
A c3 
A ci 
A c5 
A sl 

A s2 



1 



4(1 - 9a) 
3 

2(1 -5a) 
1 



1 - a ( 21 - 18A 2 + —A 4 
5 



9A 2 3 

7T 

1 - 9a 



A 



6-3 



(1 -5a) 2 
18a(l - 9a) 

(l-5a) 2 
_2a(l -9a) 2 

(l-5a) 3 

A 2 3 
4(1 -9a) 2 

A^A^^y 

2(l-5a)(l-9a) 
4 2 



(l-5a (3-2A 2 ))-^p 



(1 - 9a) 2 - A 2 mn)2 



9a 



5a) 



(B.10) 

(B.H) 
(B.12) 

(B.13) 

(B.14) 

(B.15) 

(B.16) 

(B.17) 



and the expressions of A^a and A^g.^2 are given in Eqs. (3.23)-(3.24). Note that, up to 
the overall factor 1/c 4 ,, which set the amplitude of the trispectrum, we expressed Eq. (B.9) 
in terms of c s , and not c-p, and of e (2.11). The practical reason is the following: upon 
using a redefined variable 



c 



1 — 5a 



1 — 3a 

the quadratic action (3.15) can be cast in the form 



c 



s^ 2 



dt dTX r-^ 



i2 

c 



(B.18) 



(B.19) 



which is formally equivalent to the quadratic action in terms of ( in k-inflation [41]. Results 
about the trispectrum derived in this context [35] can thus be directly translated here 
provided we use the variables ( (B.18), e (2.11) and c s (3.16) in intermediate calculations. 
Using the form 

5a 1 



S. 



(3) 



off 



, l3 a 3 eH 2 M 2 1 
dt d 3 x — 1 



o , A^n :i - A 2 c 2 ^ ( ^ )27 r^£ (B.20) 

c% 1 — 3a Cj,(l — 9aJ a 1 

of the cubic action Eq. (3.22), together with Eqs. (B.9) and (B.19), we thus straight- 
forwardly deduce from [35] the final result Eq. (3.33), except for the contact-interaction 
trispectra Eqs. (3.34)-(3.35) - coming from the higher-derivative operators 7r 2 ((did 'jit) 2 — 
(dfn) 2 ) and (dk7i) 2 ((didjii) 2 — (d 2 n) 2 ) absent in k-inflation - that we calculated from first 
principles following Eq. (B.4). 
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C Shapes of trispectra 



In this appendix, we represent the trispectrum (3.33) for the interesting value a = 0.097 
lying in the middle of the la region favored by observational constraints on the bispectrum 
(see Fig. 2), and at which the correlation of the bispectrum with the orthogonal template 
is the largest (see Fig. 1). We set cx> = 1 without losing generality and follow the same 
strategy as in [25]: we represent 2-dimensional slices of relevant trispectra in different 
representative limits, as was originally done in [35]. Note also that we chose to plot the 
scale-independent quantities T = T l{k\k 2 k 3 k 4 f'^ rather than T itself. The four limits 
we considered are: 



• The specialized planar limit, in which k\ = k 3 = ku, and the tetrahedron reduces 
to a planar quadrangle with [35] 

1/2 

(C.l) 



• 12 



k ' + ^kJ ( hh + V ~ k ^ Ak " ~ k ^) 



Shapes are then represented as functions of k 2 jk\ and k±/k\. 

• Near the double-squeezed limit: k 3 = fc 4 = k\ 2 and the tetrahedron is a planar 
quadrangle with [35] 

A /p (U2 T p _|_ p\ _ 1.2 p I p p ipp i p p _ p p _ p i pp 

■y ft<l ^ "'3 "4J "sl"s2 ' 12 14 ' "'12 4 ' "14"4 "14"3 "4 ^ "3"4 



V^2fc 4 

(C.2) 

where fc s i and k s2 are defined as 



fcs! = 2a/ (fcife4 + ki • k 4 )(kik4 - ki • k 4 ) , 

= 2 v /(fc 3 fc4 + k3-k4)(fc3fc4-k3-k4) • (C.3) 

Shapes are then represented as functions of k±jk\ and k\±/k\. 

• The folded limit: k\ 2 = 0, hence k\ = k 2 and k% = k±. Shapes are then represented 
as functions of k^/ki and k^/ki, and we assumed k& < k% without loss of generality. 

• The equilateral limit: k\ = k 2 = k% = k&. Shapes are then represented as functions 
of k\ 2 jk\ and k u /ki. 

In Figs. 4, 5, 6 and 7, we represent the (rescaled) shape function T in Eq. (3.33) for 
a = 0.097, and for comparison the representative 'equilateral-type' shape T cl currently 
used in data analysis [37], in the specialized planar limit, near the double-squeezed limit, 
in the folded limit and in the equilateral limit respectively. From these plots, the (ap- 
proximate) proportionality between the two trispectra is obvious. More generally, we have 
checked that T c \ is a good ansatz for the DBI Galileon trispectrum (3.33) in the whole 
interesting range 0.095 < a < 1/9. 
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(a) \tNL — t*Ni\l '^nl in the 3<r region allowed by the bispectrum con- 
straints, for a > 0.095. 
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(b) \tNL — £jvxl/Atjvx in the ltr region allowed by the bispectrum con- 
straints. 

Figure 3. In the plane (op, a), confidence level \tjsrL— £nl\/ ^nl at which DBI Galileon inflation 
is excluded (or not) from the observational constraint on the trispectrum. 
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Figure 4. We consider here the specialized planar limit with k\ = k^ = k±4, and plot T c \ (left) 
and T in Eq. (3.33) for a = 0.097 (right), as functions of k2/k\ and k^jk\. 




Figure 5. We look here at the shapes near the double squeezed limit: we consider the case 
where k% = k^ = kyi and the tetrahedron is a planar quadrangle. We plot T c \ (left) and T in 
Eq. (3.33) for a = 0.097 (right), as functions of k^/ki and k\±jk\. 



Figure 6. We consider here the folded limit k\2 = 0, and plot T c % (left) and T in Eq. (3.33) 
for a = 0.097 (right), as functions of k^/ki and k\^jk\. 
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